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The classical dynamics in stationary potentials that are random both in space and time is studied. 
It can be intuitively understood with the help of Chiriltov resonances that are central in the theory of 
Chaos, and explored quantitatively in the frameworl?; of the Fokker-Planck equation. In particular, 

#vj ' a simple expression for the diffusion coefficient was obtained in terms of the average power density 

,_^ , of the potential. The resulting anomalous diffusion in velocity is classified into universality classes. 

("^ ■ The general theory was applied and numerically tested for specific examples relevant for optics and 

^\| ' atom optics. 
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I. INTRODUCTION 



Potentials that are random both in space in time are models of noise. Often one seeks for methods of suppressing 

noise, the present work on the other-hand is focused on interesting properties of dynamics generated by noise. It was 

I— I subject of many sophisticated studies for nearly 100 years [l|-|4|. The response to forces resulting of such potentials 

1-^ . typically differs from ordinary diffusion. Specifically, the diffusion coefficients predicted by such mechanisms sensitively 

Sn ■ depend on the velocity of the particles. Also, if the potential is time-dependent, the energies of the particles will not be 

a constant. The existence of such 'anomalous diffusion' has been demonstrated for classical dynamics with spatially and 

temporally fluctuating potentials |5l4lC|. These works claim universal behavior in the sense that for generic random 
potentials the diffusion coefficient has a universal power-law dependence on velocity v, such that D{v) ^ \v\~^ as 
|w| — > oo. This in turn implies that asymptotically in time the average velocity satisfies (v'^) ^ t^'^. The average 



. , displacement satisfies (x^) ~ t^'^/^ for one-dimensional systems [5l-[qjl(j| (faster than ballistic) and (x^) ~ t^ (ballistic 
transport on average) for systems with dimension higher than one [6|. In the present paper we will demonstrate that 
this picture should be extended, and we will introduce such an extension. 
'T' ■ The present work is motivated by experiments in optics and in ator n. op tics. The random potentials induced by 

n^ \ optical waves in photonic lattices [ll| and in atom optics experiments [1^,113], r ely on transforming an intensity pattern 
C^ I into an effective potential for the light [li!, fl5| (the former) or the cold atoms [ij, [ij] (the latter) . Such potentials are 
O naturally described in terms of the Fourier spectrum of the waves inducing them, and their spectral coefficients are 
1^1 . assumed to be independent random variables. In the regime where the wave-packet has already acquired large values 
of velocity, one expects classical mechanics to provide a reasonable approximation. Therefore, we will explore here 
T— I \ the classical dynamics for potentials, defined by their spectral content or more precisely their average power spectral 
^ ■ density [16|. In particular, we consider the classical dynamics of a particle in potentials that are random both in space 
^f[j ' and in time, emphasizing the spreading of the velocity acquired by the particle, as time evolves. 
~^ In the present paper we will study dynamics of classical particles in potentials which are defined in terms of their 

yr^ Fourier components. Two types of potentials will be explored: 

(3 ; V^^^ (x, t)= I V (k, w) exp i (k • X - ujt) dkdw + c.c. (1) 

(N 



and 



y(2)(x,t)= v^'H^,t)\ (2) 



where V (k, w) is a random field chosen such that the distribution of V'^^'' (x, t) is stationary both in time and space. 
For time dependent potentials the velocity of the particles may grow, since energy is not conserved. The dominant 
mechanism of the growth of the velocity is via Chirikov resonances, namely, resonances between the particle dynamics 
and the external driving. Those resonances occur when the phases in ([T]) are stationary. For a given instantaneous 
velocity v = x, this happens for 

k • V = w, (3) 

and a more complicated expression is found for the potential ^. Indeed, transport takes place only in such regions in 
phase-space where the density of the resonances is non-zero. The Chirikov resonances provide an intuitive picture for 
the understanding of transport in phase-space. Since, the potentials are given in terms of their Fourier components it 
is relatively easy to calculate the spectral content or average power spectral density (PSD), S (k, oj), (see (|113p ) and 



with its help also the diffusion coefficient. The spread in velocity is then calculated using the Fokker-Planck (FP) 
approximation, assuming that the decay of the potential correlations is sufficiently rapid. The random potentials of the 
form ([1]) and ([2]) are natural in optics realizations since these are formed by a superposition of light beams. In recent 
experiments in optics [111. Il7|, light propagates paraxially in a disordered potential, which is produced by utilizing 
the photo-sensitivity of the medium |14l . il5] . In the paraxial approximation a probe beam satisfies a Schrdinger like 
equation where the refractive index plays the role of the potential, the propagation direction plays the role of time 
and the dispersion relation is approximately lo (k) = fc^/2. Since the index of refraction is proportional to intensity 
of the electric field of the writing beam, that varies like V^^^ (x, t) of (H]), the resulting potential will have the form 
of F'^-* (x, t) ([2|). Potentials of this type are also relevant for atom optics. There the electric field induced a dipole 
moment in a neutral atom which interacts with the electric field and results in a potential which is proportional to 
the intensity of the electric field. The time dependence is induced by a modulation of the intensity of the electric 
field, such that it is independent of k. 

Although this work was motivated by experiments in optics, where the dynamics is wave dynamics, the classical 
dynamics of particles in such potentials is a fundamental question by itself. Therefore, in this work we will study the 
classical dynamics of particles in potentials of the form ([T|) and ([2|). Nevertheless, since it is generally believed that 
for large velocity (or short wave length) wave dynamics is approximated by classical dynamics, we expect at least 
qualitatively to describe the original wave problem. 

We will classify the various systems into universality classes according to the velocity dependence of the diffusion 
coefficient in large and small velocity limits. The regimes of the validity of the FP approximation will be studied 
as well, and the crossover between uniform acceleration and diffusive behavior (in velocity) will be studied. The 
formalism will be demonstrated for specific examples of the form ([T|) and © . Our main result is that although the 
universality of the form D (v) ~ v^^ holds always for dimensions higher than one, for one dimensional systems there 
is a large variety of novel universality classes. A classification scheme for such classes will be presented and it will be 
demonstrated how these can be realized. 

In Sec. In] the time required for the particles to 'discover' the fiuctuations of the potential will be evaluated. After 
this time we expect the dynamics to be controlled by the FP equation. An expression for the diffusion coefficient in 
terms of the spectral content or more precisely the average power spectral density (PSD) is derived in Sec. IIIIl it is 
computed explicitly and its significance is discussed for specific examples. A scaling property of the FP equation is 
presented in Sec. IIVI A detailed calculation of the diffusion coefficient in velocity for dimensions two and higher is 
presented in Sec. IV Al for potentials of the form of ([1]), and in Sec. IVBI for potentials of the form ^. In Sec. IVII the 
rich variety of the universality classes found for one-dimensional systems is presented and the asymptotic expansion 
in powers of v~^ is discussed. The results are summarized and discussed in Sec. IVIII 

II. TIME DOMAINS 

In this section we discuss the typical time scale to observe the anomalous diffusion. For this purpose we assume 
an initial distribution that is narrow in phase-space. Initially the velocity is growing linearly in time (acceleration), 
since it takes time for the fluctuations of the potential to become effective. First we define the characteristic length 
and time scales of the potential. Consider the potential given in ([1]), let l^ be the length over which the potential has 
an almost constant derivative, resulting in a constant force. It satisfies. 



In a similar way we define the characteristic time scale Zj, 



{V')-{V)\ (4) 



It 



{v^) - {VY 



(5) 



\ ((f) 
Using, 

{V^) = / dcjidwa / dkidk2/F(ki,a;2)V'*(k2,a;2)\expi((ki - k2) • x - (wi -^2)^), 

= I Alo f dkS{k,uj), (6) 



one finds, 



W\^\ ^ / '^'^i'^^2 / dkidk2ki • k2 /V"(ki,a;2) V* (k2,iiJ2)y expi ((ki - k2) • x - (wi - ^2) t) , 

= /dw f dkk^S{k,u;), (7) 

where 5 (k, oj) is the average power spectral density (PSD) defined in (fTB]) . Therefore the characteristic length scale, 
Ix, is given by, 

,2_ /d../dk5(k,^) ^g^ 



=" /da;/dkfc2S'(k,a; 

and similarly the characteristic time scale. It, is given by, 

^ /d../dk5(k,^) 
* /dc^/dkw2^(k,u;)" ^' 

Note, that ^t and Ix are different from the exponential decay rates of correlations. The difference follows from the 
fact that we are interested in the typical lengths in both time and space where the force is constant, and not in 
the asymptotic behavior of the correlation function. For times t < It the force resulting from a potential is almost 
constant, provided that the particle is not displaced more than Ix- For a weak force, 



^0 < 7i, (10) 

the particle will travel a distance less than l^ in time It, therefore the cross-over out of the uniform acceleration regime 
takes place at time, 

U ~ If (11) 

If the force is strong so that (fTO|) is not satisfied, the particle will reach ,1^ at a time shorter than It and therefore t* 
will be, 

U - M. (12) 

V ^0 

In this work we will consider initial conditions of particles with a narrow distribution both in velocity and position, 
such that 5v < Ix/h, where 5v is the width of the distribution centered at u = 0. For these initial conditions uniform 
acceleration will take place for time t < It, while for longer time-scales the velocity will exhibit anomalous diffusion. 
In this work we will assume that the force is weak, such that (ITUl) holds and the crossover to the diffusive regime will 
take place at ([TT|). 

III. ANOMALOUS DIFFUSION IN VELOCITY 

In this section the variation of the velocity of particles by potentials V (x, t) , random both in space and time will 
be studied on time scales longer than i*. The equations of motion are, 

^ = ~vy(x,i) (13) 

dx 

d^ ="' 

where v and x are the instantaneous velocity and position. For simplicity we set the mass to be one. We will assume 
a stationary potential both in space and time and an isotropic distribution of the potential. Therefore, its correlation 
function is 

(y(xi,ii)y(x2,t2))=C(xi-X2,ti-t2), (14) 



and we may take the average of the potentials, that is a constant, to vanish, namely, 

(^(x,i))=0. (15) 

The angular brackets (.) denote the average over the ensemble of random potentials. The Weiner-Khinchin theorem 
implies, 

C (x, t) = duj dkS (k, uj) exp i (k • X - ujt) , (16) 

where S (k, uj) is the average power spectral density (PSD) of the potential, that is its spectral content. 
The Fokker-Planck (FP) equation for the velocity is given by [J|, 

where P (v, t) is the probability density and Dij is the diffusion tensor, given by the expression, 

D.,{v) = \j^ {F,{^t,t)F, {0,0)) dr, (18) 

where 

F == -Vy (19) 

is the force. Here one assumes that the correlations decay sufficiently fast and the force is sufficiently weak, so that 
the variation of the velocity on the time-scale of decay of correlations is negligible, and 

Xl-X2-v(ti-t2), (20) 

this is the usual approximation used to derive the Fokker-Planck equation [4]. Here and in the rest of the paper 
summation over repeated indexes is assumed. We will now calculate the force correlation function, 

K,j (Xi - X2,ti - ts) = {F, (Xi,ti) Fj (X2,f2)) , (21) 

which can be written using the potential correlation function, 

if.,(xi-X2,il-t2) = ^^^^^(i^(xl,il)y(x2,i2))=-^^^|x=x,-x.. (22) 

Using (frot one finds, 

Kij (xi - X2,ii - i2) = / "^^ / dkkikjS {k,uj)expi [k • (xi - X2) - uj {h - ^2)] , (23) 

and the diffusion tensor is, 

A,(v) = iy^^if.,(vT,r)dr. (24) 

Which can be written as, 

^ij i"^) = o / dr duj dk kikj S (k, oj) exp ir {k ■ V ~ uj) . (25) 

Taking the integral over r first, gives a delta function 

Aj(v) = TT I duj I dkhkjS{k,uj)S{k-v-uj). (26) 

Then the integral over uj can be easily performed to give, 

D,j{v)=n I dkkikjS(k,k-v). (27) 



Since S (k, w) is integrable by definition, we will define a function F {k,uj), so that 
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5(k,c.) = — F(k,c^). 



Then, the diffusion tensor can be written as. 



(28) 



D^J (V) 



dvidvj 



I{v). 



where 



I{v)=7T /"dkF(k,k-v), 
and we have used the fact that the integral is isotropic. Equation ([29]) can be decomposed to, 

D^J (V) 



d^ J, X ^ _d_ fv2_d]_ 
dvidvi dvi \ V dv 



U^UUj 



SijV — ViVj /v\ dl ViVjd^I 



dv v'^ dv"^ 
For a general isotropic diffusion tensor of the form, 

A,(v)=(%-5)/t(.) + ^/pW, 

the operator on the RHS of ([T7)) is. 



d 



d 



.-(d-l) 9 ^(d-l) 



d 



Lo = j-D., (v) — = „-(^-i) ^C^-D/p (v) 

ovi ovj ov ov 

We turn now to justify this result. It is convenient to use hyper-spherical variables, such that. 



d dv d 



E 



d Vj d 



dvi dvi dv ^-^ dvi dd^ v dv 



E 



d(j)j d 
dvi 



J UlfJ] 



(29) 



(30) 



(31) 



(32) 



(33) 



(34) 



where d is the dimension of space and (f)j are the angle coordinates. Using ([32|) and the new variables the operator 
reduces to, 



in - 



A 

dvi 

_d_ 
dvi 
_d_ 
dvi 
d V. 






-5f)fT{v)+^-^fp{v) 



dvi v dv 

Which using the chain rule can be written as, 

' 5iiV-ViVi/v 
^0 



\^n - -^) h («) + ^/p iv) 



[Vi - Vi) -/t {v) + —fp {v) 
V V 

d 



(v) 


d 


dvj 


(v) 


Vj d 
v dv 


d 




dv 





(35) 



d d d 

dv dv dv 

i(d-l)/p(t-)|- + |-/p(«)|- 
V ov dv ov 



Therefore, combining ((3T|) and ([33)) one identifies fp (v) = d^I/dv^. Therefore, the operator Lq is given by, 

Lo^v-^'-'^^v^'-'^Div)-^. 
dv dv 



(36) 



(37) 



The diffusion coefficient D (w) is found to be, 






D{v)^^. (38) 



It can be simplified using the definition (150)) and (pS)) . 

D(w)=7r /"dk (k-t))^5(k,k-v). (39) 



Using p9p and ahgning the a:~ component of k with velocity we get 

r'27r /'Oo 



D{v)^-^ I d9 dfc /f^+i cos^ eS {k, kv cos 9) , 
2 Jo Jo 

where Sd is the surface of the d— dimensional hyper-sphere. Changing variables to, y = w cos 9 we get 



(40) 



D{v) = ^ dy ^ y / dkk''+^S{k,ky), (41) 



where the factor of two was eliminated due the multiplicity. This expression assumes implicitly that d > 1. For one 
dimensional systems there is no integration over the angle (see Eq. (|96p in Sec. IVip . The different dependence on the 
velocity of (|4T|) and (|96p is the origin of the richness of the one dimensional behavior. For large velocities (|4T|) has 
the following asymptotic behavior, 

D{v)^^, (42) 

with 

Z?3 = 2Sd / dy / dA; y^k'^+^S (fc, fcy) . (43) 

Jo Jo 

Therefore, for dimensions higher than one the asymptotic behavior of the diffusion coefficient is indeed universal for 
any choice of the PSD and is given by (|42l) . The diffusion coefficient for zero velocity using (|40|) is given by, 



D{{)) = ^ dkk''+^S{k,0). 
4 Jo 



(44) 



If D3 y^ 0, as is the case for d > 2, the Fokker-Planck equation for the velocity (fT7|) is asymptotically equivalent to 
the equation. 



which has the scaling solution 



f - ("-"-'l'''"-'$J^ 1 -' («' 



^ <".') = 7;^9 (t I ^ <*' 



The resulting growth of the mean kinetic energy is, 

i (v^) ^ e/\ (47) 

This behavior is considered in the literature as universal [5|, |7|, |8|, |lO| . For dimensions two and higher this is indeed 
the case, however for one dimensional systems other behaviors are possible, as will be demonstrated in Sec. IVIl 



IV. SCALING PROPERTIES 

Let us assume that the average power spectral density, S (k, w), has natural spatial and temporal frequency scales, 
fcp and Wq, such that. 



V^ ~fV 



OJ 



5(k,w) = -^5(-l,-), (48) 



where Vq is a constant which determines the strength of the potential. Then ((39)) takes the form, 

D (v) = TT-^ f (k . vf S f ^, ^) dk, (49) 

fcgwo J \ko ujo J 

where 5 is a dimensionless PSD. Rescaling the variables 

k = k'/co V = v'f;o vo = -—, (50) 

ko 

gives 

D(v')^ttV;^^ /"(k'-{i)'^(k',k'-v')dk'. (51) 

Wo J 

The Fokker-Planck equation for the velocity (J17p is invariant under the transformation of variables, 

V ^- v' 

,^f[.,s^y\ (52) 

V. SPECIAL POTENTIALS 

In this section the spreading in phase-space is studied for specific potentials. Potentials of the form ([T]) will be 
studied in subsection IV Al while potentials of the form ([2]) will be studied in subsection IV Bl In this section the 
emphasis will be on the behavior for dimensions higher than one, while one dimensional systems will be analyzed in 
the next section. 

A. Potentials which are a superposition of waves 

Natural potentials to consider are potentials which are composed of a superposition of standing waves (c.f., Eq. 
©), 

T/(x,t)= /T>(k)expi(k-x-tj(k)t)dk + c.c. (53) 

with a dispersion relation uj (k) . We will assume that the amplitudes and the wave numbers are independent random 
variables. Leading to. 



where 



V (ki) V* (k^)) = Vo'f (ki) S (ki - k2) , (54) 

V (k) ) = Vq and / (k) is the distribution of the wave numbers. The correlation function of the potential 



is given by. 



C (xi - X2, ii - t2) = Vi I f (k) expi (k • (xi - X2) - w (k) {h - is)) dk + c.c. (55) 

Using (dHI the PSD is. 



S (q, Lo) — j—^ / C (x, t) exp — i (q • x — tot) dxdi (56) 

= TTTd+T^o' / / (k) exp * ((k - q) • X - {lo (k) - oj) t) dkdxdt 
(ztt) J 

+ jAd+i^o j f (k) exp -z ((k + q) • X - (u; (k) + uj) t) dkdxdi. 

By making the integrals over x and t we get, 

Siq,Lu) = Vo^fiq)S{LU-L,iq)) + Vo^fi-q)S{LU + L,iq)). (57) 

Starting from (|4T|) and assuming that both / (k) and uj (k) are isotropic gives, 

^ („) ^ flTL / dj/ ^ ^ / dfc k^+^f (k) [S [ky - LO (fc)) +5{ky + Lo (fc))] . (58) 



We can now perform the integral over y and obtain 



where il (.) is a Heaviside step function, which follows from the restriction of the delta function. In the limit of large 
velocities we have, 

D{v)^^, (60) 

with 

/■oo 

Di = 2SdV^ dkk''-^f{k)uj''{k). (61) 

Jo 



For zero velocity using (|44)) we have. 



D{Q) = ^v;^ r dkk''+^ 



f{k)5{u{k)). (62) 



- 2^°yo k(o)r^^-^' 

where we have assumed that 

LO (fc) = 0, (63) 

has only one solution fc = 0, which is a general property of physically relevant dispersion relations. 

We will proceed by calculating the diffusion coefficient for a particular choice of the dispersion relation, u (k) — fc^/2, 
and a uniform density of wave numbers, 

{-rj^ < k < kjf 
n ^ . . ' (64) 

fc > fcfl 

where Vd is the volume of the d— dimensional unit hyper sphere. Therefore, we have 

^T/2 f fn" dfc , '''''^^ V < kD/2 

^("^-2fc^.3 f^-dfc , ^""^ «>fcH/2' ^^^^ 



> 

V 



O 




-10^ 



log^o(t) 



Figure 1: A log- log plot of average squared velocity as a function of time for a two dimensional system with the potential 
(|53|l and a distribution of wave- numbers (|64p . The blue dots represent the result of the Monte-Carlo solution of (|13|l averaged 
over 20 realizations and the black solid line is the numerical solution of the Fokker-Planck equation for the velocity (|17[1 . The 
dashed black and red lines are guides for the eye with the corresponding slopes of 2 and 2/5. The initial condition was a narrow 
distribution of velocities for the Fokker-Planck and x = v = for the Monte-Carlo calculation. The parameters used for this 
simulation are, Vb ~ 10~^, kn = 0.1. 



where we have used the fact Sd/Vd = d. Taking the integrals gives, 



Div) 



_dvl_ 



d+3 12 127 2 ' 2 ' ItP 



V < fcfl/2 

V > kR/2 



where 1F2 {x) is Gauss's hypergeometric function. For the special case of d 



D[v) 




(!?) 2-^3 



\ll) 



1 



V '2v / 



(66) 



V < kR/2 

V > kR/2 



(67) 



The prediction of the FP equation was compared to the Monte-Carlo simulation and is presented in Fig. [T] It is 
found that uniform acceleration takes place for times v < kR/2 and the FP prediction clearly fail in this regime. This 



can be expected since from (p7)) D (v) 



the assumption that the velocity is constant (1^ cannot hold, rendering 



the FP approximation inconsistent. Nevertheless, for v > kR/2 the FP predictions are satisfied. Note that, for any 
potential discussed in this section the scaling properties of Sec. llVl hold. 

Numerically, it is found that for large velocities (v > kR/2) the average square position is growing ballisticly, 
(x^) ~ t^, as expected from previous studies [a,HHll3i while for small velocities (v < kR/2), (a;^) ~ i^, as expected 
for uniform acceleration. 



B. Potentials proportional to the intensity of a superposition of waves 



In this subsection the diffusion coefficient for potentials of the form ([2]) will be calculated. These potentials appear 
in experiments with neutral atoms and in some experiments in optics, where the refractive index (which plays the 
role of the potential) of a photosensitive material is proportional to the intensity of light [ij, llH. The complex field, 
[/, which represents the superposition pattern of waves is given by. 



/ 



f/(x,i)= / dkt/(k)expi(k-x-w(fc)i) 



(68) 



10 



where tj (k) are random Fourier coefficients, and w (fc) is some dispersion relation. Following the analysis of the 
previous subsection we will assume that U is composed of plane-waves with independent phases and amplitudes, 
leading to 

(69) 



C/(k)^ - 

C/(ki)C/*(k2)) =/o /(ki)<5(ki-k2) 

;7(ki){7(k2)\ = 0, 



where 



U (k) ) = Iq and / (k) is the probability density of the wave numbers k. The potential is proportional to 

(70) 



the intensity of the U, 

V (x, t) - \U (x, t)f = /"dkidks U (ki) U* (ka) expi [(ki - ka) • x - (c^ {k^) - u (fcz)) t] , 
Using the assumptions ([69]) we get that the potential is constant on average, 

(F(x,i)) = /'dkidk2 (?7(ki)C/*(k2))expz[(ki-k2)-x-(c^(fci)-w(A:2))i] 
= /o|dk/(k)=/o. 
The correlation function of the potential is, 

(y(xi,<i)y(x2,i2)) - /"dkidk2dk3dk4 (f7(ki);7*(k2);7(k3)C/*(k4) 

X expi [(ki - k2) • xi - {(jJ {ki) - to (fc2))ti] 
X exp i [(ka - k4) • X2 - (w (/ca) - w (fc4)) 12] ■ 



(71) 



(72) 



Since we have assumed that the complex amplitudes of the plane-waves are independent random variables we can 
decompose. 



C/(ki)C/*(k2)C/(k3)t/*(k4)^ = ^t/(ki)C/*(k2)^^C/(k3)t/*(k4)^ 

+ (?7(ki)C/*(k4))(f7*(k2)f7(k3)' 

= Ilf (ki) / (k3) 5 (ki - k2) <5 (k3 - k4) 
+ Ilf (ki) / (k2) 5 (ki - k4) <5 (k2 - k3) . 

Therefore, the correlation function reduces to, 
C(xi -X2,ii -^2) = ^0 

+ /2y"dkidk2/(ki)/(k2) 

X expi[(ki -k2) • (xi - X2) - {uj {k{) ~ uj {k2)) {h - ^2)] • 



(73) 
(74) 



The PSD is just. 



S{ci,oj) = - — -^^ / dx / di C (x, t) exp -i (q • X - wi) 



(27r) 
/2<5(q)<5(c.) 



(75) 



A_ |dxy" di y dkidk2/ (ki) / (k2) 



(27r) 
X expi [(ki — k2 — q) • X — [w (fci) — cj (^2) — cj] t] . 

Taking the integral with respect to x and t we get, 

S (q, w) = 1^6 (q) 6 {u) + II j dkidka/ (ki) / (k2) <5 (ki - k2 - q) <5 [w - (u; (fci) - w (fe))] , 



(76) 

(77) 
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^(q,c.) ^ I^S icO S {u;) + 1^1 dkfik)fik -(i)S[io- {CO ik)-uji\k- q\))]. (78) 

The diffusion cocfBcient is given by (see, ([39l) ) 

D (v) = TTii J dqdk (q • vf / (k) / (k - q) 5 (q • v - (c. (fc) - c. (|k - q|))) . (79) 

Changing the variables to, 

q' = k q 

k' = k, (80) 

and suppressing the primes for simplicity gives, 

D (v) = 7r/2 I dqdk ((k - q) • vf f (k) / (q) 5 [(k - q) • v - (c. (fc) - c. (<?))] . (81) 

Choosing the x- components of k and q such that they are ahgned with the velocity v yields, 



2 /iQo r2TT />oo r2-K 

2^. 



D{v)^^ll[^] / dqj dej dfc / d9k{qkY'-' (82) 



'0 "'0 ^0 ^0 

%) f (k) f {q) 5 (W (fc Cwo Vfc y wo ^q 

Changing variables to, yk = v cos 9k and yq — v cos 9q in (|5^ we have the Jacobian, 



X {k cos 9k — q cos 6q) f {k) f {q)S {v {k cos 9k ~ q cos 9 q) — {uj{k) — u; {q))) 



J= ^ ^. (83) 



v^^[l-{yk/vf){l-{y,/vf) 



In these variables the diffusion coefficient is, 



D iv) = 47r/^ ( ^ ) 4 / dg / dfc / dy^ / dy, (qfc)""^ (84) 



'S'd\ 1 / j_ / _,,_ / j_ / j_ ^_,_Mi-l 



Jo 

2 



y(i-W«)')(i-(y»') 

For large velocities this gives, 



^^^' '^^'^ f{k)fiq)S{{kyk-qyq)-{ujik)-u;{q))) 



D{v)^^, (85) 



with 

Ds = -I'.SJ I dq q''-^ I dfc fc"-^ (c. (fc) -uj{q)rf{k)f (q) . (86) 



A /'OO /'OO 







Since / (fc) is non-negative the integral for D^ does not vanish and for an appropriate choice of / (fc) it is convergent. 
The diffusion coefficient for zero velocity can be calculated from (|82| by setting v = 0, 



/ r< \ ^ /"CxD /'27r /"Oo /'Zir 

i^(0)=<f^jy dqj d9qj dkj d9k{qkt~^ (87) 



X {kcos9k - q cos 9qf f{k)f{q)S{u}{k) - cj (q)) . 



Assuming that lo (fc) is a monotonic function we have. 



S{cj{k)~u:{q))^j^^S{k-q), (88) 



12 




log^o(t) 

Figure 2: Same as Fig. [T]but with the potential (|70p and a uniform distribution of wave- numbers given by (|64p . The dashed 
black lines are guides for the eye with the corresponding slopes of 2 and 2/5. The parameters used for this simulation are, 
Jo = 10-^ kR = 0.1. 



and therefore, 



^(9)1' 



D{0) = ^Io(^) J d0, /" d^fe (COS0, 

/•OO 

X / dq (J 
After integrating over the angles, 

Jo \^ Wl 



k — cos ( 



(89) 



(90) 



that is non- vanishing. Therefore, the potentials described in this section are generic, namely, produce the universal 
behavior both for small and large velocities. We will now demonstrate the calculation for a particular choice of 
dispersion relation, uj (k) = k'^/2 and a uniform distribution function / (fc) on a disk given by (|64p the integrals in 
(1821) cannot be taken explicitly, however we have calculated them numerically (see Fig. [SJ and the resulting dynamics 
is demonstrated in Fig. [51 Also here we found numerically that (x^) ^ t^ for large t and (x^) ~ t"^ for small t. 
For same dispersion relation and a Gaussian distribution of wave numbers. 



using (|8T|) we have, 
D{v) = 



nil 



{2nklf 27r 



/(k) = 



i2nklf' 



dr / dqdk ((k— q) ■ ij) e ^'°« expir 



(k- 



where we have used the exponential representation of the delta function. Noticing that. 



D{v) 



nP, 



{2TTk' 



2tt 



dr 

'72 



dqdk e '^''r — — ^ exp it 



(k 



{k'-q') 



{k'-q') 



(91) 



(92) 



(93) 
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and taking the Gaussian integrals gives 

D{v) 



{2TTklf 27r 



°° dr (2^fc^) 



(1 + k%T^f'^ dv- 



fc2 ^2 

-e \i+%^'''' 



(94) 



i?(«)=^/2 



dr 



271 J_^ T^ [I + k%T^f^^ dv"^ 



(95) 



The dynamics fohowing from (j95p is rather similar to that demonstrated in Fig. [2] and therefore it is not presented 
here. 



VI. ANOMALOUS DIFFUSION IN ONE DIMENSIONAL SYSTEMS 

In this section we will consider one dimensional systems and explore the different possibilities of asymptotic behavior 
of the velocity distribution. In particular, we will show that the asyniptotic behavior of the form (j42l) . is not universal 
since some possibilities were not considered in previous studies [5|, |8|, |lO| . In the first subsection we will present a 
some representative examples of the new universality classes, and in the second subsection we will explain when the 
asymptotic expansion of [5|, |8|, llO| is not effective, and provide a prescription for designing transport properties using 
the PSD. 



Following from 



A. New universality classes 

the diffusion coefficient for a one dimensional system takes the form. 



D{v)=n f k'^S{k,kv)dk. 



For a PSD with the property, 

S (fc, Lo) — uj/k > v„ 
implies that the diffusion coefficient is zero for large velocities, 

D{v) ^0 V > Dmax 



(96) 



(97) 



(98) 



as demonstrated in [lali^l- The correlation function corresponding to the PSD with the property ([W)) may be infinitely 
differentiable, nevertheless the diffusion coefficient vanishes for large velocities as is clear from (|98)) . a possibility which 
was overlooked in 0, |^, [T^] . 

We will now present specific examples of this new universality class using the specific potentials of the form ^ and 
dU. The PSD of potentials of the type ^ is. 



S {k, u) = V^f {k) [6{uj-LJ {k)) + 6 {lj + oj (k))] . 
Taking the dispersion relation uj — fc^/2, and using ([M]) and ([M|) yields, 

D (v) = nVi I k^f (fc) [5 {kv - k^/2) + 5 {kv + k^/2)] dfc 

"(5(fc) S(k-2v) 5{k + 2v) 



(99) 



nV,' / k'fik) 



dfc, 



(100) 



which after integrating over k gives. 



Div)=87TVo^\v\fi2v) 



(101) 
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Any asymptotic behavior is possible as the diffusion coefficient is proportional to / {v) and therefore clearly may decay 
faster than v~^ . The correlation function for these potentials can be obtained for example for a Gaussian distribution 
of wave- numbers ([9T|). Using ([T6| and (|99|) we have, 



cos 



Cix,t)^2Vo^ 






2(l + fei*^) 



i arg (l + ik%t) 



{l + k%t^f' 



exp 



2{l + k%t^) 



(102) 



Note, that this correlation function is an infinitely differentiable function over the whole {x, t) plane. It has fast 
decaying spatial correlations compared to a very slow temporal correlations. 
For potentials of the type (O the PSD is given by, 

(103) 



S {k, u) = /o^ j dqf (q) f{q-k)S[iJ-iLJ {q) ~ UJ [q - k))] . 
As in the preceding example, taking a quadratic dispersion relation, lo = fc^/2, gives 

D{v) = „ll j k'l(,)f{g-k)S t« - i (,» - (5 - t)') Aqik 



\v-q\ 



4vr/2 / \q^v\f{q)f{2v-q)dq. 



\v-q\ 



Changing the variables to q = q' + v and suppressing the prime we have, 

D {v) = 4< f \q\ f{q + v)fiq-v) dq. 



And the decay of the diffusion coefficient with the velocity is also dictated by the decay of f (q 
may be obtained for example for a uniform distribution of wave-numbers (|64|1 . giving |16l.llE 



D(v) 



'^{kn-\v\f \v\<kn 







\v\ > kR 



(104) 

(105) 
An explicit expression 

(106) 



The resulting dynamics is demonstrated in Fig. [3l The regime of a unit slope corresponds to regular diffusion in 
velocity and the asymptotic regime corresponds to an absence of diffusion in velocity. We have verified numerically 
that the asymptotic regime corresponds to a ballistic growth in position, (x^) ~ t^ ■ 

Another explicit expression may be obtained for a Gaussian distribution of wave-numbers (|91|) . The diffusion 
coefficient is given by. 



D{v) = 2/^e ^. 



(107) 



The predictions of this equation are compared with Monte-Carlo simulation and the results are presented in Fig. 3] 
The correlation function for this potentials can be obtained using (fTBI) . 



C(x,t)^ll 



yTTfcp 



exp 



1,2 2 

1 + 4*2 



(108) 



Similarly, to the correlation function (I102p this correlation function is infinitely differentiable, decaying fast in position 
and slowly decaying in time. 

This existence of a dispersion relation between the spatial and temporal frequencies is not a necessary condition for 
the new universality class, as is demonstrated by the following PSD, 



S{k,uj) 



eVS 



47rfcowo 
where the normalization factor was chosen such that 



■e V »= / e 



4- + -^ 

^0 ^'^o 



/ 



S (/c, w) dfcdw = V^ 



(109) 



(110) 
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log^o(t) 



Figure 3: A log-log plot of average squared velocity as a function of time for a one dimensional system with a potential (|70|) and 
a distribution of wave-numbers ((MJ. The blue dots represent the result of the Monte-Carlo solution of p3p averaged over 20 
realizations and the black solid line is the numerical solution of the Fokker-Planck equation for the velocity (|17|) . The dashed 
black and red lines are guides for the eye with the corresponding slopes of 2 , 1 and 0. The initial condition was a narrow 
distribution of velocities for the Fokker-Planck and x = v — for the Monte-Carlo calculation. The parameters used for this 
simulation are, Jo — 10~*, kn = 0.1. 




-10 r- 



Figure 4: Same as Fig. |3]but with a distribution of wave- numbers (|91[1 . The dashed black lines are guides for the eye with the 
corresponding slopes of 2 and 1. The parameters used for this simulation are, Vo ~ W~ , kR = 0.1. 



Using (j96p and defining vq = cjo/fco yields. 



Div) 



eVJ fcp 

\/7r Vq 



-^Jl + {v/vo)'^ 



(111) 



which is exponentially decaying for a large velocity. We could not obtain an analytical expression for the correlation 
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function, but it can be expressed as an integral, 

C(x,t) = ,° / e-iTJ e V^ ^>J expi(kx-ujt)dkdio, 
AnkoujQ J 

The temporal correlations are rapidly decaying and since the function in the integrand is infinitely differentiable and 
integrable following from the Riemann-Lebesgue lemma the spatial correlations will also decay faster than any power 
law. 

In summary, in this subsection we have presented a few examples of potentials, which give diffusion coefficients 
with various asymptotic behavior for large velocities different from (|^^ . We have also calculated their corresponding 
correlation functions and showed that they are smooth and differentiable. Therefore, this class of potentials spans 
a new universality class, which is qualitatively different from the universal behavior described in previous studies 
p, [a, |l3| • We have also demonstrated that the conditions for appearance of this class are not necessarily connected 
to long range correlations or the special form of the potentials ([T]) and ^. In particular, the existence of a dispersion 
relation is not required. In the next subsection we will show that the conclusions drawn here are valid for any physical 
dispersion relation, and present a classification of the potentials, based on their PSDs, into different universality 
classes. 

B. Classification of potentials into different universality classes 

In the present work we have chosen to work with the PSD and not the correlation function, since in relevant 
experiments [111 . IitI . I19| the PSD rather then the correlation function is naturally controlled. For the purpose of 
clarity we will repeat here the asymptotic expansion of the diffusion coefficient, which was done in [8|, |lO| and connect 
it to our representation using the PSD. 

The PSD is the Fourier transform of the correlation function of the potential, namely, 

S{k,uj) = 2 / C {x,t)exp—i{kx — u}t)dxdt. (113) 



D(v) ^ — / fc^C (x, t) exp -ik (x - vt) dxdtdk, (114) 

47r 



Therefore using (|96 



which could be written as, 

1 r d^ 

D(v) = C(x,t)— ^ exp -ik(x-vt) dxdtdk. (115) 

An J ox'' 

Integrating by parts twice and taking into the account that C (x, t) and its derivatives vanish in the limit of a; — > ±oo, 
we get, 

1 f d^C 

D (v) — / (x,t]exp—ik(x — vt) dxdtdk. (116) 

47r / ox'' 



Finally integrating over k gives, 



D{v)^-^ f ^{x,t)5{x-vt)dxdt. (117) 



This is the same expression for the diffusion coefficient, which was obtained in previous studies (see Eq. 7 of ^^5)- For 
simplicity of notation we will set. 
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Changing variables to t — t' /v we get, 



D{v) = -^ fw(x,-]6{x-t')dxdt' 



Therefore, if W is sufficiently differentiable, the asymptotic expansion of D (v) is. 



Div) = ^ + ^ + ^---, (120) 



W V v^ 



where 



1 f 1 d'^'^W 
Dn = -^x^'-'-g^;^{x,T)\r=odx n=l,2,...,. (121) 



The first term vanishes, since 



Di = -\jw{x,{))dx = -]^j^{x,Q)dx = 0, 



(122) 



and all the even terms vanish due to the symmetry, x — >■ — a; (following from the fact that the potential is translationaly 
invariant). 

Therefore, the leading term in the asymptotic expansion of the diffusion coefficient is, 

D{v)^^, (124) 

where 

D3^--J x'^-^ix,r)\r^odx. (125) 

In [3, 1^; [l3 it ^^s claimed that if the correlation function is at least twice differentiable in both arguments then the 
asymptotic behavior (|124p follows, and furthermore it is universal. However, as was demonstrated in the previous 
subsection exceptions are possible, since some or even all of the terms I?„>3 may vanish, giving a qualitatively different 
asymptotic behavior. If the diffusion coefficient vanishes or decreases faster than any power law for large velocities 
the asymptotic expansion in powers of v~^ clearly vanishes. Thus the transport in such case is sub-diffusive for large 
velocity and the growth of the moments of the velocity distribution is slower than any power law. In what follows we 
will demonstrate that for the potentials of type ([T]) and ([2]) the asymptotic expansion of D (v) vanishes, for any choice 
of the dispersion relation, provided that the distribution of the wave-numbers decays faster than any power law. 
Taking (|M1) and the relation (I113P we can calculate, 

W{x,t) = -V^ f k^f{k)[S{uj-uj{k))+5{uj + uj{k))]expi{kx-ujt)dkduj 

= -Vi^ k'^f{k)[expi{kx~uj{k)t) + expi{kx + uj{k)t)]dk. (126) 



Using the expression for the coefficients of the asymptotic expansion (J12ip we have, 

D„ ^ Y^ [ ef{k){-iuj{k)Y'^t"'-^e'''''dkdt' 



n 
2 r 



1/2 r 8^-^ 



ef [k) c.«-i (fc) jr—^e^'' dkdt'. (127) 



18 

Integrating by parts n — 1 times, and assuming that / (fc) decreases faster than any power law, so that the boundary 
terms may be disregarded, we obtain. 



n! / dk 



Dn = ' ', ° e^'^* -^(fcV(fc)-"-Mfc))dfcdi' 



2«-l , .2 



27r 



{-lY^'-'V^ 



d 



n-l 



9A:' 



— (fc2/(fc)a."-i(fc)) 



fc=0 



For a physical dispersion relation which has the property uj {Q) — Q we will have I?„ = for all n. As was shown in 
(jlOip the diffusion coefficient in this case is not zero, but a function which decays with velocity faster than any power 
law, rendering the expansion (|120p useless. For a potential with a PSD given by ()103|) we have, 

W (x, t) = -/q / fc^/ [q) f {k - q)S[uj - (w (g) ~uj{q ~ k))] exp i (/ex - tut) dqdkduj 

= -/q / k^f {q) f{k-q) exp i {kx - (cj (q) - w (q - k)) t) dqdk. (128) 



The coefficients of the asymptotic expansion are therefore given by, 

Dn = ^jH{u:{q)-Uj{q-k))r-'k'f{q)f{k-q)t'"-^e^'^''dt'dqdk 



^ y (-» (^ (<Z) - c^ (9 - A:)))"-' fcV (?) / (fc - 9) (^-*^ j e^'* di'dgdfc. 



(129) 



Integrating by parts and assuming / (fc) decays faster than any power law gives, 






e 



° 2d y aA^ 

2n! aA^ 



n— 



- {[c. (g) - c. (9 - fc)]" ' fcV (g) / (A; - g)} df'dgdfc 



2vr/g ^ ,1, ^,„,, ^ / [c.(9)-c.((7-fc)r'^V(9)/(fc-?)d4 . (130) 



fc=0 



These coefficients vanish for any dispersion relation, and the diffusion coefficient is different from zero, and was 
obtained explicitly for a specific choice of the dispersion relation in (|105|) . 

In experiments one usually does not control the diffusion coefficient directly, and therefore we will classify the 
potentials leading to slow transport using the PSD, which can be controlled. Starting from a relation for the diffusion 
coefficient (|96l) we change the variable to fc = k' jv and obtain. 



D{v) = ^ fk'^s('^,k'\dk'. 



(131) 



If the PSD is not singular (note, that this is not the case for the potentials of type ([TJ and ([2]) discussed earlier) and 
differentiable on the line A: = (where k is the first argument of S (A;, a;)) we can expand, 

D{v)^^ + ^ + ---, (132) 

V V 



where 



{n-Sy. J " dk' 



Dn ^ ,^ ^„ / dc^c^"-^^^^^^ {k,Lo) Ifc^o n>3. (133) 



Note, that in this representation the leading order of the expansion is directly v~^ , since representation using the 
PSD already contains the assumption that the potential is stationary. If the first nmax derivatives, d"S/dk", vanish 
on the line A; = 0, where rimax is the first derivative that is different from zero, then the resulting asymptotic behavior 
is D {v) ^ Dn^^.^+s/v'^""^^^^. All the derivatives vanish only if S{k,uj) is non-analytic on the line A; = 0. The 
non-analytic behavior may result either from a function, which is strictly zero on some finite strip around the line 
A: = 0, as in ([W)) or due to an essential singularity of S{k,uj) on this line, which will be demonstrated bellow. The 
first case will lead to a diffusion coefficient that will vanish for large velocities, for example poip (with finite support 
/ (A:)) or (|106p . The second case will result in a decay faster than any power law, for example (J107I) . In experiments 
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with a good control of the PSD, by controUing the amount of vanishing derivatives of the PSD on the hne k — 0, one 



can vary the asymptotic behavior of the velocity distribution. The range of variation is from D (v) 



, through 



D (v) ~ v~" (with n > 3), sub-exponential, exponential, super-exponential and up-to D (v) = (for v > Umax)- 

We will now turn to examine the correlation function of the potential, C (x, t), for potentials, which belong to the 
new universality classes. For simplicity we will consider the behavior along the line. 



c(i) = C(0,t) 



S {k,uj) dfc ) e 



and, 



'duj 



S{k,uj)dLu]e'^''dk 



c{x) = C{x,0) -- 
The moments of the correlation function along these lines are, 

(r) = frc{t)dt 
and 

For the PSD 1^ we have, 

(i") = V^ ff'f{k)[S{u;-u;{k)) + S{uj + ujik))], 



(134) 



(135) 



(136) 



(137) 



Mfcdwdi 



Vo I f (k) t" (e-^"('='* + e*"''^)*) dfcdi. 



(138) 



We will now calculate the second moment, which will contain all the representative features. Using the fact that, 

32 



{t')^2Vo'lfik)(^-^yosu;t 



dfcdi. 



and 



we have. 



Integrating by parts twice gives. 



it') 





d dk d 






duj doj dk ' 




'^^"7 ^^''\du:dkd^dk) 


cos wt dfcdi. 


r/adfca fdk ^ 
^^" J [dkdujdk[dJ^''>, 


1 J cos wt dfcdi 


471^0^ 


If d dk d /^dfc ^ 
\ dk duj dk \dio j 


0] 


fc=0 



(139) 



(140) 



(141) 



(142) 



where we have used the assumption that a; (fc) = if and only if fc = 0. Note, that for any increasing (near fc = 0) 
dispersion relation, except a; = |fc|c, this moment diverges. For the later dispersion relation the derivative, dfc/dcj 
does not exist at fc = 0. Therefore, the temporal correlation function c (i) for these potentials is slowly decreasing. 
Nevertheless, for the spatial correlation function, we have. 



(143) 



h(k) ^ / S'(fc,w)da; 
= 2yoV(fc)- 
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Figure 5: Log-Log plot of the diffusion coefficients as a function of tiie velocity for a two dimensional system and different 
potentials. The dots stand for potentials of the form ([l]), while the solid lines stand for potentials of the form ((2]). The blue 
(light) line and dots are for Gaussian distribution of wave-numbers while the black line and dots are for uniform distribution 
of wave-numbers over a disc of unit radius. The dashed black lines are guides to the eye with the slopes, 2 and —3. 

If one chooses the function / (fc) to be infinitely differentiable then following from the Riemann-Lebesgue lemma the 
spatial correlation function, c (a;) will decay faster than any power law. A similar calculation could be repeated for 
the potentials of the type ([2]), with the PSD (jl03p . leading to a similar result. Therefore, a natural question to ask is 
whether a necessary condition for the new universality class is a slowly decaying correlation function. The answer to 
this question is negative as was exemplified after (|109p . 

In this subsection we have demonstrated that potentials of the type ([T]) and ([2]), with any dispersion relation and 
a fast decaying distribution of wave- numbers will result in a transport slower than (jl24[) . Those potentials will have 
short correlations in position and long correlations in time. For potentials with smooth PSDs we have suggested a 
classification scheme based on their PSD properties. 



VII. DISCUSSION 



In the present work the diffusion coefficient for the velocity is presented in terms of the spectral content, namely, 
the average power spectral density (PSD), (see Eq. (p9)) ) for potentials that are random both in space and time. This 
representation is very natural for potentials generated as a superposition of waves used in optics and atom optics. It 
is of particular use in the case where the potential is a stationary process both in time and space. The simplicity 
of the expression enabled to explore the properties of the diffusion coefficient, to establish a scaling form of the 
Fokker-Planck equation, and to discover new universality classes. In particular, we were able to calculate explicitly 
the diffusion coefficient for representative examples, relevant for applications both in optics and atom optics (Gaussian 
and uniform distributions of wave- vectors). We have shown that for dimensions larger than one only one universality 
class is possible in the framework of the FP approximation, D (w) ^ w~'^. However for one dimensional systems new 
possibilities for large velocity asymptotics were also found. 

In the limit of small velocities the diffusion coefficient vanishes only when S (fc, 0) ~ 5 (/c), which is not the typical 
behavior. However, for potentials described in Sec. IV Al Z? (0) = and moreover for small velocity the diffusion 
coefficient grows faster than v, invalidating the Fokker-Planck approximation for those potentials, in this regime. 

There is an important case were the waves composing the potential satisfy a dispersion relation given by to (fc) ~ k^ . 
In this case we found that the diffusion coefficient depends on the characteristic scale of the potential, kj^. In Fig. [5] 
D {v) is presented for similar values of the characteristic scale of the potential, kn- 

The main result of the paper is the classification of universality classes in one dimensional systems presented in 
Sec. IVIl In the past it was found that in the large velocity limit the diffusion coefficient depends on the velocity as, 
D (u) ~ v~^ [5|, |6|, |8|, [lOJ. In the present work we have shown that this is always the case for dimensions larger than 
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one, and explained the mechanism of this behavior (see (|4T|) ). However for one dimensional systems this is only one 
of the possibilities. Generally it depends on the asymptotic expansion of D [v) in powers of v~^: 

1. The first term of the asymptotic expansion of D (v) is non-zero, D {v) ^ v~^ . 

2. The first non- vanishing term in the asymptotic expansion of D (v) is the n— th term then, D (v) ^ v~^. 

3. All terms in the asymptotic expansion of D (v) are zero, D (v) < v~", for any a > 0. In particular, the diffusion 
coefficient may be zero for v > Wmax (e.g. (|106p ) or non-zero but decreasing faster than any power law (e.g. 

dMD). 

All these possibilities can be realized in experiments with control over the PSD or more precisely / (k) , as is clear 
from poip or (llOSp . Unlike previous studies 0, S Hy] the classification into universality classes does not rely on 
the differentiability of the correlation function of the potential. Additionally, as was demonstrated in (|112[) . the new 
universality classes does not depend on the range of the correlation function. The general behavior can be understood 
heuristically in terms of the Chirikov resonances [20| . 

We have also shown that initially particles will experience a uniform acceleration and at some later stage, on a 
longer time-scale, the velocity will follow diffusion as predicted by the Fokker-Planck equation. Eventually, it will 
reach the asymptotic long-time behavior predicted by the Fokker-Planck approximation. 

In the present research the spreading of the velocity distribution is studied in the framework of the Fokker-Planck 
equation, therefore an obvious question to study is what happens when the Fokker-Planck approximation fails. An- 
other important issue, which was not addressed in this work, is how to obtain analytically the spreading in position. 
Here this spreading was studied numerically for the potentials (H]) and ^. Since in the experiments that have mo- 
tivated this work, the relevant dynamics is of waves, rather than particles, an obvious question to explore is the 
correspondence between the classical and wave dynamics. 

This work was motivated by the experimental work of Liad Levi and Mordechai Segev, whom we thank for many 
stimulating discussions and for providing crucial insight for this problem. It is our great pleasure to thank Tom 
Spencer for introducing us to 0, 0, [l3| and to Michael Wilkinson for introducing us to [3, [l3| • Many of the results 
of the present work originated from fruitful discussions with Michael Wilkinson during his visit to the Technion. We 
would also like to thank Igor Aleiner, Boris Altshuler and Michael Berry for informative discussions. This work was 
partly supported by the US- Israel Binational Science Foundation (BSF) , by the Minerva Center of Nonlinear Physics 
of Complex Systems, by the Shlomo Kaplansky academic chair, by the Fund for promotion of research at the Technion. 
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